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Based on the dissipative Landau-Lifshitz equation, the spa- 
tiotemporal structure formation problem is investigated in the 
region far above the transverse ferromagnetic resonance insta- 
bility. Apart from the external fields, the model contains an 
isotropic exchange field, a shape demagnetisation field and an 
anisotropy field. Numerical simulations exhibit in the rotating 
frame a stationary domain structure with a precessing motion 
in the walls regimes. Employing analytical methods, charac- 
teristic elements of this structure are explained. This driven 
dissipative system shows similarities to equilibrium systems 
of coexisting phases and organises itself in such a way that 
the local dynamics tends to become hamiltonian. 
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Ferromagnetic systems strongly driven by external 
transverse magnetic fields have been under investigation 
for more than four decades (for the early work see 
Based on a nonlinear spin-wave expansion around the 
homogeneous ground state, Suhl Q showed that cer- 
tain modes become unstable for increasing amplitude of 
the pump field. Detailed measurements exhibit beyond 
threshold many of the phenomena predicted by the gen- 
eral theory of non-linear dynamics. Based on ||^, trun- 
cated spin-wave mode models have been published, which 
successfully describe the experimental findings (for a re- 
cent, comprehensive presentation see [^). 

Apart from single exceptions like the simulations of 
Elmer Q|, who found dynamic domains for a special 
model, the existing investigations are limited to the weak 
nonlinear regime. The regime far above the instabilities 
with pump amplitudes of the order of the demagnetisa- 
tion field remains an open question (compare |^ p. 1075). 

This question is of high interest in the general theory 
of spatiotemporal pattern formation in dissipative sys- 
tems [||. At the beginnings of this modern and inter- 
disciplinary field of research Anderson |^ proposed that 
a driven ferromagnet is the characteristic example for 
pattern formation and for the general problem of how 
the concepts of equilibrium phase transitions can be ex- 
tended to driven dissipative systems. Recently, Cross and 
Hohenberg Q raised doubts against this proposal. 

It is mainly the dipole interaction which makes pattern 
formation in ferromagnetic materials a difficult problem. 
The difficulties already arise in the undriven case . In 



a first step, the theory of static ferromagnetic structures 
Q usually employs an approximation. Above all, stray 
fields are omitted and only the part of the dipolc interac- 
tion which describes the shape demagnetisation fields is 
considered. Even in this approximation many interesting 
results have been obtained Q both for the formation of 
static structures and for the wall dynamics. 

Using this approximation for the dipole field, this work 
adresses the pattern formation problem in driven mag- 
netic systems. Focusing on the region far above thresh- 
old, it is the aim of the present investigations to work out 
general features and mechanisms by employing both ana- 
lytical and numerical methods. It is beyond the scope of 
this paper to improve the existing near threshold treat- 
ments, as this would certainly require the complete dipole 
field used in the spin- wave approaches. 

At a mesoscopic scale the dynamics is governed by the 
Landau-Lifshitz equation as was recently demonstrated 
by microscopic investigations pHlO|]. In the frame rotat- 
ing with the driving frequency oj around the direction, 
this equation of motion takes the form 

m = -mx(H'=ff -cjG;,) -Fm X (m X H°^) . (1) 

m(r, t) and H"^^ are the local magnetisation and the effec- 
tive field in the rotating frame, being related to the quan- 
tities miab and Hf^ in the laboratory frame by miab = 
exp(wte2x) m and by Hj^^ = exp(wte2X ) H°^, re- 
spectively. F represents the Landau-Lifshitz damping 
rate. In reduced units the gyromagnetic ratio and the 
magnitude of the magnetisation m = |m| are equal to 1. 

Specifying the model under investigation, it is assumed 
for the effective field that 

H'^* = h\\ez + /^_Le^ - m + JAm + Am^Gz (2) 

where and are the amplitudes of the external static 
and the external circular driving rf field, respectively. 
The term m = V^^ J m dV represents, again in reduced 
units, the demagnetisation field of a sphere of volume V. 
The contribution JAm results from the isotropic ferro- 
magnetic exchange interaction. An uniaxial anisotropy 
is described by Anizez- 

In the first part of this letter the results of numerical 
simulations in one spatial dimension with an arbitary ^ 
direction m(r, i) — > m(^,t) are reported using periodic 
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boundary conditions m(^, = m.{^+L,t). These investi- 
gations were performed on a vectorized super computer. 
The program uses an Euler integration scheme in time 
and is based on semi-spectral methods. 

These simulations exhibit a temporal evolution of 
m(^, t) which is characteristic for structure formation 
in dissipative systems. After a period of transient be- 
haviour, which strongly depends on the special initial 
state and which may be very complex, a tendency to- 
ward the formation of domains is found. In the early 
stage these domains interact which each other, merging 
into domains of larger size. With increasing time and in- 
creasing domain sizes, this process slows down and finally 
the structure becomes stationary on large time scales. 

Fig.l shows such a final domain state obtained by a 
long-time study starting from a randomly disturbed, ho- 
mogeneous initial state. The magnetisation is nearly ev- 
erywhere constant, taking the values m+ or m_, respec- 
tively. Contrary to this behaviour within the domains, 
a strong time dependency is found in the narrow wall 
regimes which separate the domains. As Fig. 2 and Fig. 3 
demonstrate, the magnetisation m at a fixed position is 
anharmonically oscillating with a characteristic internal 
period of rotation Tpci-, which is found to be indepen- 
dent of the specific position. This spatiotemporal wall 
structure is dynamically very stable. In the simulations 
no changes of any significance could be found for times 
corresponding to 10^ periods of Tpor. 

It should be pointed out that this wall structure as 
well as the macroscopic quantities m+ , m_ , m and Tpor 
are independent of the specific initial state and thus are 
characteristic elements for the structure formation. 

To analyse these characteristic features found in the 
simulations a perturbative approach has been worked out 
using the multiple-time scaling method and treating the 
contributions J Am and Am^e^ as perturbations [prl . 
The theory has been worked out up to the first order. 
The first order treatment is lengthy and rather technical 
thus I restrict myself here to the zeroth order. 

In this order eqs.(|^) and (||) reduce to 



m(r,t) = -mxHi - Fm X (m X H ) . (3) 



with 



H2 = Hi +0162 = h\\ez + h±ex - m(t). (4) 

Due to the rn term, the problem described by eqs.(||) and 
is of the mean field type and the usual technique can 
be applied. For the present case this would imply solv- 
ing eqs.(H) for given rnit) as a first step. Integration of 
the obtained solution m(r, t) over the sample volume V 
would then lead to a self-consistency condition for rn{t) , 
from which rn,{t) and the complete solution can in prin- 
ciple be obtained. This procedure can in general not be 
performed, as already the analytic solution of the first 



step is not known The fixed points of eqs.(|^), how- 
ever, and their stability can explicitly be determined by 
applying the mean field procedure. 

The results of this, basically straightforward, fixed 
point analysis show that marginally stable solutions ex- 
hibiting a domain state structure are possible. The lo- 
cal magnetisation m(r) takes only two values m+ and 
m_ realised in the generally disconnected partial vol- 
umes V+ = uj^V and V- — n^V, respectively. 

In such domain states the fields Hi and H2 become 
orthogonal 
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and it is advantageous to write the results in terms of the 
internal coordinates introduced by ei = Hi/_ffi,e2 = 
H2/H2 and by 63 — 61X62- The transformation from 
the internal to the primary frame is found to be 



(6) 



where 

d^hii-io; M=(d2+r2)i/2 (7) 
v=iu^-r^hl)'/^ ; w = iu' + hiy^' (8) 

were introduced. The magnitudes of the fields Hi and 
H2 are calculated to take the values 

Hi = ujuw^^ ; H2 — u)h±'w~^ (9) 

and the domain magnetisations are found to be given by 

m± = u-^Tvei + Thj_e3). (10) 

For the bulk magnetisation in a domain state, the self- 
consistency condition leads to 
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and the n± are calculated to be 
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The domain states exist only in a subspace of the param- 
eter space spanned by , h±, lo and by F. This subspace 
is implicitly determined by the condition < n+ < 1. 
For parameter values outside this subspace eq.(|^) has a 
different stable fixed point solution characterised by a ho- 
mogeneous magnetisation m(r) = mhom- These homoge- 
neous states, being of minor interest for the present work, 
are unstable in the regime where the domain states exist. 
Apart from the domain and the homogeneous states, no 
other stable fixed points of eq.(||) exist. 

Next, the findings of the linear stability analysis about 
the domain state fixed points are presented. Employing 
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the usual exp(A<) ansatz for the deviations, one finds for 
all modes with the exception of four 



X = ±i fl with Q 
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which implies that they are undamped and oscillat- 
ing with f2. The four remaining modes are of collec- 
tive character as they describe disturbances of the form 
Jy^ (m(r, t) — ni±) d^. The characteristic equation 
for these collective modes is found to be given by 



- 2n+n^V%\u^^ [5ft(A+A_) - n+nJ\ , (14) 
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where A± = (A =F ifl)/(r + i) + n± was introduced. It 
has been shown analytically for the weak damping case 
(r ^ 1) and numerically for the general case, that the 
collective modes are damped everywhere. Note that this 
implies relaxation of rn{t) to mdom ■ 

Under the constraint rn{t) = rndom, the analysis of 
eqs.(|^) and (jj) can be extended to the nonlinear regime. 
Employing eq.(H), this leads to 
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ujw ^ (^u mxei + rh±^ m x (m x e^)^ . (15) 

For this equation of motion the quantity 

TOi(r,t)(r/i_Lm3(r,i) - u)-^ = vC{r) (16) 

is a constant of motion and cq.(p^ is integrable, as in 
addition mf + + m; 



2 + mg = 1 holds. Eq.(|16[) represents a 
plane and thus the resulting motion of m(r, t) is a preces- 
sion on a cone at an arbitrary but fixed position r. This 
precession is found to be anharmonic with the period of 
rotation of Tpcr — 2TTil~^ independent of r. 

The values of C(r) are restricted by the condition 
|C(r)| < 1, as otherwise the plane ([l^ ) does not intersect 
with the sphere. For the limiting values C(r) = ±1 the 
plane becomes tangential and consequently m becomes 
time independent taking the fixed point values m± of 
eq. (0). Thus a wall is characterised by a change from 
C=— ltoC=l. Considering now the ^ oo limit, 
the constraint rn{t) = mdom can be satisfied, as long as 
the entire wall volume goes to zero faster than V. 

Summing up the analysis of eqs.(^ and (^ has shown 
that stable solutions of domain type do exist. In the do- 
main regions the local magnetisation is stationary but 
precesses with Tpcr in the wall regime. It is obvious that 
these results are qualitatively in accord with the numeri- 
cal findings. Using the parameter value of the simulation, 
quantitative agreement is found for mdom, for m±, for 
n± and for Tpor- The deviations are less than one percent. 
Recall also that the numerical results of Figs. 1-3 were al- 
ready presented in terms of the internal coordinates (|^). 
Furthermore, the numerical results show that the motion 



in the wall regime is planar and well described by eq.(16) 
Finally, the first order treatment ||l^ of eqs.(|l|) and (|^ 



is able to explain with high accuracy the complete spa- 
tiotemporal wall structure of Figs. 2 and 3. 

Some features of the present investigation which may 
be of general interest should be pointed out. 

The behaviour of this dissipative system driven away 
from equilibrium in the stationary domain state is analo- 
gous to a thermodynamic system at a first order equilib- 
rium transition exhibiting coexisting phases. The trans- 
lational symmetry is broken for the domain structure. 
Macroscopic variables, as in the present case mdom, rn± 
and n±, are independent of the initial conditions and 
only depend on the external parameters like h± and 
u!. In this context eq.(||) has to be interpreted as the 
criterium for existence of the domain states. Note that 
this criterium reduces for the static case (w = 0) to 
/lyez -|- h±e,j: ~ rrldom = 0, which means that the in- 
ternal field vanishes. As this is the usual criterium for 
a ferromagnet to be at the equilibrium phase transition, 
some aspects of the present work can be interpreted as an 
extension of the concepts of phase transitions to driven 
systems. 

Another interesting aspect of the present work is the 
local generally nonlinear dynamics within the domain 
regimes. As long as the sample is sufficiently large, this 
local dynamics is integrable or is hamiltonian according 
to the investigations of eq. (^5|) . Although the system be- 
comes weakly dissipative in the first order treatment , 
it is worthwhile to point out this tendency of the system 
to organise itself on the global macroscopic scale in such 
a way that the local mesoscopic dynamics becomes at 
least nearly hamiltonian. 

In this context it is interesting to go back to the mi- 
croscopic scale where the system is a quantum me- 
chanical many body spin system in contact with a bath 
of sufficiently low temperature. It is basically the compe- 
tition between the short ranged exchange interaction and 
the long ranged repulsive dipole interaction, which leads 
to the formation of structure and to the nearly hamilto- 
nian dynamics on the mesoscopic scale. 

Ferromagnets are representative of a whole class of 
many body systems characterised by short ranged attrac- 
tive interactions in competition with long ranged repul- 
sive interactions. Thus, the question arises whether the 
tendency toward a local hamiltonian dynamics is generic 
for all these rather common systems. Should this be the 
case, it would perhaps become clear why structure and 
why hamiltonian dynamics are so widespread in nature. 

Coming back to the ferromagnet, it is pointed out that 
neither the location nor the geometric shape of the walls 
is determined by the present approach. By analogy with 
the static case it is expected that the dipolar stray 
field and the usual boundary effects will reduce this free- 
dom. The formation of regular patterns in driven ferro- 
magnets seems to be possible. 

For completeness, the results of this work are spe- 
cialised to the weakly nonlinear regime <C 1. For 
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simplicity, the discussion will be restricted to w = h\\ 
and to r ^ 1. A critical value of the pump field 
l^cnt ^ Y{hj - 1)1/2 results. For h± < h"'* the homoge- 
neous state is stable and for h± > h"'^*- the domain states 
arise. As in this work the dipolar stray field is neglected, 
^crit (jjffej-s from the threshold values of |^ , which under 
usual conditions are smaller than h'^"^. Lastly, the find- 
ings for the power absorption P will be given. For the 
homogeneous state P = T~^h']_ and for the domain states 
P = r/i^/i^(r2 + /i^)^i are obtained. This result implies 
an experimentally well established [|l]j3| saturation of P 
for h± ^ r. Thus, as a side product of this work it must 
be concluded that this saturation is not solely explained 
through spin- wave approaches. 

I am particularly indebted to F. Matthaus for his pro- 
gram from which the numerical results have been ob- 
tained. I gratefully acknowledge discussions with H. Bon- 
ner, H. J. Elmer, W. Just and G. Sauermann. This work 
was performed within SFB 185. 



FIG. 1. Spatial dependence of a stationary domain state 
in the rotating frame. Plotted are cosi^ = m\{rn\ + m^)~^^^ 
and ma, where rrii are the components of the local magneti- 
sation in the internal coordinates, defined by eq.(6). The 
simulation uses 1024 mesh points and the parameter values 
= = 2, ftx = .5, J = .01, A = -.005, r = .1 and 
L = 1024. 

FIG. 2. Spatiotemporal structure of cos cj> in the vicinity 
of a wall located at ^waii with Tper = 18.6. Definitions and 
parameter values as in Fig.l. 

FIG. 3. Spatiotemporal structure of mz in the vicinity of a 
wall. Definitions and parameter values as in Figs.l and 2. 
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